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We write a Renormalization Group (RG) equation for the function / in a theory of gravity in the 
f{R) truncation. Our equation differs from previous ones due to the exponential parametrization 
of the quantum fluctuations and to the choice of gauge. The cutoff procedure depends on three 
free parameters, and we find that there exist discrete special choices of parameters for which the 
flow equation has fixed points where f — fo -\- fiR -I f 2 R^■ For other values of the parameters the 
solution seems to be continuously deformed. 
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Introduction. One way to achieve a quantum- 
field-theoretic UV completion of General Relativ¬ 
ity is to find a fixed point in a “theory space” 
parametrized by the most general diffeomorphism- 
invariant semilocal functional of the metric. The 
study of the RG flow in this infinite-dimensional 
space is a daunting task. A convenient simplified 
environment where one can study the RG flow of 
gravity is given by effective actions of the form 

Tk = j d^x^fiR), (1) 

where / depends on the RG scale k. This is called 
the f{R) truncation. The first studies focused on 
the case when / is a polynomial [1-4], and a fixed 
point was found there. This study has now been 
brought to high order [5]. A RG equation for / 
had been written early on, but serious attempts 
to find functional fixed points (henceforth called 
“scaling solutions”) were not made until [6]. Neg¬ 
ative results concerning the equations written in 
[2, 4] have been reported in [7, 8]. It has been 
argued that the solution may require truncations 
that go beyond functionals of the background field 
alone, and taking into account the Ward identi¬ 
ties of the quantum-background split symmetry. 
This is a general issue that goes beyond the /(i?) 
truncations and progress in this direction has been 
made in [9, 10]. In the meantime solutions were 
found in simplified (lower-dimensional and/or con¬ 


formally reduced) settings [11, 12] and recently also 
in the full four-dimensional case [13]. 

Here we use in the f{R) truncation the lessons 
learned in scalar-tensor theories with arbitrary 
scalar potential V{(j}) and nonminimal interactions 
F{(j))R. It was found there that the use of the 
exponential parametrization for the metric and of 
a “physical gauge”, based on removing the spin- 
zero and spin-one gauge degrees of freedom of the 
quantum fluctuation, simplifies the equations to 
the point where finding exact scaling solutions in 
closed form becomes possible [14, 15]. While we do 
not expect the classical equivalence between f{R) 
and scalar-tensor theories to extend to the quan¬ 
tum domain, these results encourage us to try a 
similar strategy also in the case of the f{R) trun¬ 
cation. 

Flow equation for f{R) Gravity. We use the 
background field method. Motivated by the non¬ 
linear nature of the space of metrics we split [16]: 

: ( 2 ) 

where g is the background field and h the quantum 
field. Henceforth we assume that the background 
space is a sphere. The York decomposition = 

4“ ^ + V -\- V^V,ytT — (T -|- -^g^ii/h, 

where = Oj leads to the 

following very neat form for the off-shell Hessian: 
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Here s = h — V^cr is the gauge-invariant combina¬ 
tion of the spin-zero degrees of freedom. The only 
gauge-non-invariant degree of freedom appearing 
in the Hessian is the trace part h, and it multiplies 
the classical equations of motion. 

Since we are interested in the off-shell beta func¬ 
tions, this structure suggests removing h as a gauge 
choice. This then leaves a structure that is the 
same on- and off-shell. Actually, if we work on 
the sphere, there is not enough gauge freedom to 
remove the constant mode of h, which remains as 
a single, isolated quantum mechanical mode. One 
important virtue of the exponential parametriza- 
tion is that results are independent of the way 
in which one fixes the residual volume-preserving 
diffeomorphisms (with a reparametrization of the 
trace, this is true for the full diffeomorphism group 
[17]). We will discuss this point elsewhere. Here, as 
in [14], we simply choose the physical gauge = 0. 


This leads to a real transverse vector and scalar 
ghost, with Lagrangian 

The latter cancels the contribution of the same op¬ 
erator appearing in the s-s Hessian. 

To write the flow equation we employ for the 
spin-i field a cutoff of the form f{R)Rk{Ai), where 

Rf^{z) =_(fc2 _ _y2 _ 

Ai = —— yi?, A 2 = —— aR. The param¬ 
eters a, /3, 7 , introduced in [11], correspond to 
different choices of cutoffs. According to standard 
lore, if there is a scaling solution, its coordinates 
will generally depend on such choices but overall 
qualitative features should not. By the standard 
procedure, we get the functional renormalization 
group equation 


rfc= -Tr 
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where the dot denotes k-^ and Pk{z) = z + Rk{z). 
Here we have not coarse-grained the isolated mode 
of h. We will comment on this below. The traces 
have to be performed summing over the spectra of 
the operators A^ using the explicit expression for 
the eigenvalues Ak and multiplicities Dn ■ The 
primes on the traces indicate that the lowest eigen¬ 
values of Ai and the two lowest eigenvalues of Aq 
are left out, so that all sums over the spectral 
values start from n = 2 [4]. The sums can be 
evaluated using the Euler-Maclaurin formula. In 
four dimensions, defining the dimensionless vari¬ 
able (p{r) = k~‘^f{R), r = R/k^, one arrives at the 
following flow equation; 
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Coarse-graining also the isolated mode h = const 

with the cutoff would produce an additional 

1 ^ 

contribution - ?■ 

1277*^ 16-\-2(p—r(p' 

The same equation can be derived using the heat 
kernel expansion. (Due to the form of the cutoff, 
only the coefficients up to 65 are needed.) We note 
that the equation is simpler than the ones that had 
been written previously. In particular, compared 
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to [13], its denominators are of lower order and 
the function ip does not appear undifferentiated 
in the r.h.s.. In the Einstein-Hilbert truncation, 
this reduces to the statement that the cosmological 
constant does not appear in the flow of Newton’s 
constant [14], a fact that allows nonsingular flows 
in the infrared. 

Scaling solutions. The normal form of the flow 
equation has a singularity at r = 0 and further 
fixed singularities depending on /3. The analysis 
of [7] showed that isolated solutions are expected 
to occur when the number of fixed singularities 
matches the order of the equation. This is the 
case when /? < 0.3945 (counting only singularities 
for positive r). Instead of analyzing numerically 
the equations for fixed a, /3, 7 , we treat these pa¬ 
rameters as unknowns to solve for. For large r the 
scaling solutions are expected to grow like , pos¬ 
sibly with logarithmic corrections. The simplest 
possible solutions are therefore of the form 

= go+ giT + g 2 r‘^ . ( 6 ) 

If we insert this ansatz in the flow equation and 
write the fixed point equation as a single fraction, 
its numerator is a fifth order polynomial in r. By 
equating to zero the coefficients of all powers of 
r one obtains a system of six equations for the six 
unknowns go, gi, g 2 , a, l3 and 7 . This system has a 
number of solutions whose properties are reported 
in the table. 


10®a 

10^/3 

10®7 

lO^go* 

10 ^ 51 * 

10^ g2* 
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-593 

-73.5 

-177 

7.28 

-8.42 

1.71 

3.78 

-616 

-70.7 

-154 

7.42 

-8.64 

1.74 

3.75 

-564 

-80.3 

-168 

6.82 

-8.77 

1.83 

3.70 

-543 

-87.4 

-126 

6.31 

-9.47 

2.06 

3.43 

-420 

-100.5 

-3.19 

4.90 

-10.2 

2.83 

2.93 

-173 

-2.98 

244 

4.53 

-8.34 

2.70 

2.18 

-146 

-64973 

250 

2.90 

-10.7 

0.0006 

2.58 

-109 

-22267 

307 

2.90 

-10.4 

0.0045 

2.45 

109 

-3564 

526 

2.84 

-7.83 

0.094 

C 

377 

-1305 

794 

2.57 

-4.37 

0.214 

> 4 


TABLE I: The properties of the exact quadratic so¬ 
lutions: parameters (first three columns), couplings 
(nest three columns) and critical exponent (last col¬ 
umn). ”C” stands for ’’complex. 

Since the denominator of the equation has at 
least one zero for positive r that is not matched 
by a corresponding zero of the numerator, the so¬ 
lutions shown are valid on the whole positive real 
axis except for isolated points. 


All fixed points have a critical exponent that is 
exactly equal to 4, corresponding to the volume de¬ 
gree of freedom. The last column reports the most 
relevant exponent besides this, obtained from an 
analysis based on a polynomial expansion around 
the origin. The first six lines are solutions with ex¬ 
actly two relevant directions; the two subsequent 
ones have three relevant directions. The last two 
lines have four relevant directions, of which two 
are complex conjugate (denoted by C); further¬ 
more the last has a critical exponent that is larger 
than four. The other solutions have only real criti¬ 
cal exponents, contrary to previous analyses in full 
f{R) gravity. We consider this a desirable feature. 

We observe that for all these scaling solutions, 
the equation of motion 2/ — Rf = 0 has the so¬ 
lution i?* = —2gok^/gi, and therefore avoid the 
possible issue of the redundancy of all eigenpertur- 
bations [ 8 ]. Since go > 0 and 51 < 0 the solution 
is compatible with spherical topology. 

The numerical similarities between the six so¬ 
lutions of the first group and the two solutions 
of the second group suggests that perhaps they 
are the same scaling solution in different cutoff 
schemes. (The quantitatively non-negligible dif¬ 
ferences by factors of order one are typical for the 
scheme dependence in this type of calculations.) 
To substantiate this hypothesis we have studied 
the polynomial approximation to the flow equa¬ 
tion along straight paths in the space of param¬ 
eters joining the different solutions. (For exam¬ 
ple, a{t) = -0.593 - 0.0227t, /3(t) = -0.0735 -f 
0.0028It,7(t) = —0.177 -|- 0.0226t interpolates be¬ 
tween the first two solutions when 0 < t < 1.) We 
have considered truncations to sixth order polyno¬ 
mials, and find that for any pair of solutions in the 
first group there exists a continuous interpolation. 
Likewise, the two solutions of the second group are 
continuously related. 

In addition to the expansion around r = 0 we 
have also considered the asymptotic behavior of 
the solutions for large r. For generic values of a, 
j3, 7 it is of the form 

(p{r) = Ar"^ + air -I- aq -I- a_i/r -|- a_ 2 /r^ + ... 

where A is a free parameter and Oi are given func¬ 
tions of a, P, 7 . For 0 = 1/6 also the parameter A 
is fixed. For /3 = 0 or 7 = 1/4 the leading term goes 
like log r and also subsequent terms have loga¬ 
rithmic corrections. It will clearly be important to 
match these behaviors and establish numerically 
the existence of scaling solutions of (4) for other 
cutoff choices. 

While we do not have an exact, nor a complete 
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numerical solution except at the special points 
listed, the existence of continuous polynomial de¬ 
formations suggests very strongly that the exis¬ 
tence of the scaling solutions is not limited to the 
special cutoff choices listed in the table. Instead, 
we conjecture that the scaling solution exists in a 
whole open subspace of parameters and that the 
effect of choosing special values of the parameters 
is to give the scaling solution a very simple func¬ 
tional form. 

We have also looked for polynomial solutions 
with more traditional cutoff choices. With the so- 
called type I cutoff a = /3 = 7 = 0, truncating at 
order eight polynomial, we find the solution 

(p{r) = 0.00482-0.00463r-H0.00131r2 
-1.86 X 10“V3 - 5.02 X 10“®r^ - 9.11 x 10"V 
-2.00 X 10“ V - 4.23 X 10“®r^ - 9.51 x 10“®r® 

with critical exponents 4, 2.21, —2.51, —5.21, 
-7.50, -9.52, -11.45, -12.9, -15.1. This can be 
seen to be a continuous deformation of the first 
group of fixed points in the table. 

In [4] a type-I cutoff was also used, so it is use¬ 
ful to compare this solution to the one given in ta¬ 
bles III and IV there. The differences, which can 
be ascribed to the different parametrization of the 
quantum field and to the different gauge choice, 
are significant. While go is almost the same, gi 
is here smaller by a factor five, which means that 
G* is larger by a comparable factor. This is in ac¬ 
cordance with the results of [14]. Whereas in [4] 
the coefficients gi remained of comparable size as 
functions of i, here they decrease quite strongly. 
This is similar to the findings of [18] in unimodu- 
lar gravity. Most significantly, there are here only 
two relevant directions instead of three, as found 
so far in all treatments of the f{R) truncation. 

In the same truncation but with the type-II cut¬ 
off a = —1/6, P = 1/3, 7 = 1/4, deforming the 
first group of solutions leads to the fixed point 

ip{r) = 0.00426- 0.00702r-h 0.00280r2 
-2.60 X 10“V - 1.31 X 10“®r^ - 7.06 x 10“V 
-h5.47 X 10“®r® -h 1.43 x 10“^r^ -b 8.31 x 10“^°r® 

with critical exponents 4, 1.91, —1.54, —4.73, 
-7.52, -10.0, -12.4, -14.3, -23.3. 

From a preliminary analysis of polynomial ex¬ 
pansions, it would seem that the solutions of the 
second group cannot be continued to type I or type 
II cutoffs. 

Alternative equation. Performing the spectral 
sums as in [6], leads again to equation (4) but with 


different coefficients: 
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This equation admits again some exact quadratic 
solutions. The interesting ones are 


lO^a 

10^/3 

1037 

lO^go* 

lO^gi, 

10^32* 

e 

-97.8 

38.9 

319 

4.31 

-7.46 

2.85 

2.03 

-438 

-122 

-21.0 

4.67 

-10.4 

3.14 

3.2 

134 

-2765 

551 

2.82 

-7.70 

0.13 

C 

505 

-715 

922 

2.16 

-2.65 

0.21 

>4 

-564 

-63.8 

-147 

7.83 

-6.80 

1.35 

>4 


The first has two relevant directions and is the 
analogue of the first six solutions of (5); the fol¬ 
lowing one has three relevant directions; the third 
has four relevant directions and the last two have 
at least one critical exponent larger than four, pos¬ 
sibly a complex conjugate pair, but with very slow 
convergence with the order of the polynomial. 

Much of what was said for the other equation 
can be repeated with small changes in this case. 

Discussion. The main positive result of this 
work is the existence of special quadratic scaling 
solutions for particular cutoff choices, which re¬ 
main close to this shape when the cutoff is var¬ 
ied generically. These solutions provide further 
evidence for asymptotic safety of gravity. The 
fact that they show similarities to the Starobin- 
sky model of inflation (which is favored by current 
cosmological data) is also encouraging. 

While these results definitely constitute progress 
relative to polynomial truncations, the situation is 
still not satisfactory. The existence of physically 
undesirable solutions with critical exponents larger 
than four suggest that these may be artifacts of the 
approximations made and that they may go away 
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in an improved treatment. More generally, the ex¬ 
istence of mnltiple solutions makes the identifica¬ 
tion of the physically relevant one not obvious. In 
this respect there is of course the ever-present need 
to go beyond the single-field truncations and to 
take iirto account the split-symmetry Ward iden¬ 
tities. But even within the single-field truncation 
there may be room for better understanding. One 
aspect that remains somewhat unclear is the phys¬ 
ical irreaning of the limit r —> oo, <C R) i.e. 
of coarse-graining at leirgth scales that are much 
larger than the diameter of the manifold. In fact 
it had been suggested in [12] that one may not 
even need to solve any equation beyond some fi¬ 
nite value of r. Here, as in [13], we do solve an 
equation for all 0 < r < oo, but for both equations 
we have considered, their form for large r is some¬ 
what ambiguous. This is because both the use of 
the Euler-Maclaurin formula and the method of [6] 
are two arbitrary ways of analytically interpolating 
the spectral sums. The existence of isolated singu¬ 
lar points in our solutions may be a manifestation 
of the same problem. We observe that these issues 
could be avoided using a non-compact background 
manifold. 
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